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Houston,  Director,  CERC. 
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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

acres 

13,277.1 

square  meters 

degrees  (angle) 

0.0174S329 

radian 

feet 

0.3048 

maters 

mile  (US  Statute) 

1.6093 

kilometers 

nautical  mile 

1.852 

kilometers 

pounds 

2.2046 

kilograms 

Aooesslon  For 
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1  Introduction 


Background 


At  the  request  of  the  U.S.  Army  Engineo'  Division,  Pacific  Ocean  (POD), 
a  numerical  model  wave  response  study  of  proposed  improvements  to  Maalaea 
small  boat  harbor  was  conducted  by  the  U.S.  Army  Engineer  Waterways 
Experiment  Station’s  (WES’s)  Coastal  Engineering  Research  Center  (CERC). 
The  study  was  conducted  to  d^ermine  an  optimal  design  plan  of  improvonent 
which  would  provide  the  harbor  with  adequate  protection  from  the  incident 
wave  climate.  The  ocisting  harbor  facility  consists  of  an  entrance  diaimel, 
turning  basin,  two  protective  breakwiUers,  93  berths,  a  haulout  and  launching 
ramp,  and  a  100,000-lb‘  crq)acity  cold  storage  plant.  Following  evaluation  of 
the  existing  harbor,  five  propos^  design  plans  of  improvement  were 
investigated. 


Study  Location 


Maalaea  small  boat  harbor  is  located  on  the  soudiwest  coast  of  the  island  of 
Maui,  HI,  the  second-largest  island  in  the  Hawaiian  chain.  The  harbor  is 
!q)proxiinately  7  miles  south  of  the  County  seat  in  Wailuku  and  approximately 
8  miles  south  of  the  commercial  and  business  center  of  Kahului  figure  1). 

The  Sordine  of  Maalaea  Bay  is  part  of  an  isthmus  connecting  two  inactive 
volcanos  which  form  west  and  east  Maui.  The  shoreline  is  characterized  by  a 
long  narrow  coral-sand  beach  and  includes  the  world-renowned  Maalaea 
Pipeline  surfing  area.  Maalaea  Harbor  is  located  at  the  extreme  west  end  of 
this  beach. 

Maalaea  Harbor  was  first  devdoped  by  the  Territory  of  Hawaii  in  1952. 
The  harbor  was  modified  in  1955,  1959,  and  to  the  present  configuration  in 
1979  (Figure  2).  The  existing  facility  consists  of  a  90-fl-wide,  12-ft-deq;> 


*  A  table  of  facton  for  convoting  non-SI  uniU  of  meaaurement  to  SI  uniU  is  presented  on 
pagev. 
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Figure  1 .  Study  location 
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Figure  2.  Existing  plan 


entrance  channel;  an  11.3  acre  dredged  basin;  a  100-ft-long,  90-ft-wide  break¬ 
water  and  revetment  on  the  south  side  of  the  basin;  an  870-ft-long  breakwater 
on  the  east  side  of  the  basin;  and  a  300-ft-long,  50-ft-wide  paved  wharf  on  the 
north  side  of  the  basin.  The  capacity  of  the  harbor  is  limited  to  93  small  craft 
due  to  hazardous  conditions.  A  100,000-lb-capacity  cold  storage  plant  was 
constructed  for  use  by  commercial  fishermen  operating  out  of  Maalaea 
Harbor.  A  launching  ramp  is  located  at  the  west  end  of  the  harbor  basin.  A 
U.S.  Coast  Guard  cutter  is  stationed  at  Maalaea  Harbor. 

Maalaea  Harbor  experiences  problems  which  include  severe  harbor  surge, 
entrance  channel  navigation  difficulties,  and  inadequate  harbor  facilities.  The 
surge  results  from  the  existing  configuration  and  alignment  of  the  harbor 
entrance,  which  allow  direct  wave  propagation  through  the  channel  opening. 
Surge  problems  cause  navigational  hazards  and  prevent  safe  berthing  in  some 
portions  of  the  harbor. 

Proposed  improvements  to  Maalaea  harbor  are  limited  by  several  con¬ 
straints.  The  most  signiftca''t  is  that  the  harbor  site  is  fixed  and  alternate  sites 
can  not  be  considered.  Constraints  mandated  by  harbor  users  and  local 
surfers  include;  (a)  the  existing  breakwater  structures  must  remain  intact  and 
changes  to  the  structures  must  be  additive;  (b)  construction  of  modifications 
must  be  accomplished  without  serious  interruption  of  harbor  navigation;  and 
(c)  additional  structures  must  not  extend  beyond  the  present  eastern  harbor 
boundary  in  order  to  avoid  impacts  on  the  surfing  area  outside  the  harbor.  In 
early  1980,  a  hydraulic  model  study  was  conducted  at  WES  to  investigate  the 
stability  of  various  breakwater  cross  sections  considered  in  the  proposed  plans 
of  improvement.  Details  of  the  study  are  provided  in  Carver  and  Markle 
(1981).  The  General  Design  Memorandum  (GDM)  for  Maalaea  Harbor  for 
Light-Draft  Vessels  (U.S.  Army  Engineer  District,  Honolulu  1980)  contains  a 
record  of  the  research  and  planning  which  led  to  proposed  design 
improvements,  Plan  1  (Figure  3). 

Plan  1  will  provide  berthing  facilities  for  approximately  310  small  craft, 
and  includes  the  following  improvements: 

a.  A  620-ft-long  extension  to  the  existing  south  breakwater. 

b.  An  additional  400-ft-long  revetment  on  the  seaward  side  of  the  existing 
south  breakwater. 

c.  A  610-ft-long  entrance  channei  varying  in  width  from  150  to  180  ft, 
and  varying  in  depth  from  12  to  IS  ft. 

d.  A  1.7-acre,  12-ft-deep  turning  basin. 

e.  Removal  of  80  ft  from  the  existing  east  breakwater  head. 

/  A  50-ft-wide,  720-ft-long  interior  revetment  adjacent  to  the  existing  east 
breakwater. 
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iyure  3.  Proposed  Plan  1 


g.  An  8-ft-deep  berthing  area  adjacent  to  the  existing  east  breakwater. 

h.  A  570-ft-long  interior  revetment  varying  in  width  from  50  to  170  ft. 

Plan  2  (Figure  4)  redirects  the  entrance  channel  to  the  west  and  includes 
the  following  improvements; 

a  Removal  of  300  ft  from  the  existing  south  breakwater  tip. 

b.  A  610-ft-long  15-ft-deep  entrance  channel,  varying  in  width  from  150  to 
200  ft,  and  varying  in  depth  ftom  12  to  15  ft. 

c.  A  1.7-acre,  12-ft-deep  turning  basin. 

d.  Removal  of  80  ft  from  the  existing  east  breakwater  head. 

e.  A  600-ft-iong  extension  to  the  existing  east  breakwater. 

f.  A  50-ft-wide,  6(X)-ft-long  interior  revetment  adjacent  to  the  existing  east 
breakwater. 

g.  An  8-ft-deep  berthing  area  adjacent  to  the  existing  east  breakwater. 

h.  A  570-ft-long  interior  revetment  varying  in  width  from  50  to  170  ft. 

Plan  3  (Figure  5)  includes  the  same  improvements  as  Plan  2  with  the 
exception  of  an  additional  extension  to  the  existing  east  breakwater.  The 
6(X)-ft-long  extension  will  continue  an  additional  250  ft  towards  the  west. 

Two  additional  plans  were  considered  in  the  second  phase  of  this  study 
(Figure  6).  P..ey  are  modifications  of  Plan  I.  Plan  la  is  the  same  as  Plan  1 
except  the  new  south  breakwater  extension  and  entrance  channel  are  rotated 
clockwise  7  deg.  Plan  lb  is  identical  to  Plan  la  except  a  vertical  sheet-pile 
bulkhead  replaces  the  revetment  along  the  east  side  of  the  center  mole. 

Headquarters,  U.S.  Army  Corps  of  Engineers  (HQUSACE)  and  POD 
established  the  following  study  objectives:  Verify  that  the  proposed  harbor 
design  improvements  meet  the  criteria  that  wave  heights  not  exceed  1  ft  in 
berthing  areas  and  2  ft  in  the  entrance  and  access  channels  and  turning  basin 
more  than  approximately  10  percent  of  the  time  per  year.  Develop  a  final 
design  plan  that  satisfies  the  locals  to  the  harbor  and  adjacent  areas.  To 
accomplish  these  objectives,  the  HARBD  numerical  harbor  wave  response 
model  (Chen  and  Houston  1987)  developed  at  CERC  was  used  to  test  the 
existing  harbor  configuration  and  proposed  Plans  1,  2.  and  3.  The  existing 
configuration  was  tested  to  establish  harbor  response  to  waves  for  existing 
conditions.  The  three  proposed  plans  were  tested  since  the  final  design  plan 
must  meet  the  aforementioned  constraints.  Subsequently.  POD  requested 
testing  of  the  two  modifications  to  Plan  1  (Plans  la  and  lb).  The  proposed 
plans  will  hereafter  be  referred  to  as  "Plan  1",  Plan  2",  "Plan  3",  "Plan  la", 
and  "Plan  lb." 


6 


Chapter  1  Introduction 


Chaptar  1  Introduction 


Figure  4.  Proposed  Plan  2 


Figure  5.  Proposed  Plan  3 


Modeling  Approach 


Both  numerical  and  physical  modeling  were  considered  for  this  study. 
Physical  modeling  has  the  advantages  of  providing  more  complete,  reliable 
results  for  this  particular  study  and  would  allow  more  comprehensive 
optimization  of  the  project  design.  However,  since  the  proposed  design  plans 
are  considerably  different,  the  physical  model  would  cost  significantly  more 
and  take  longer  to  complete  than  the  numerical  model  approach.  The  assump¬ 
tions  inherent  in  the  numerical  modeling  approach  are  as  follows: 

No  wave  transmission  or  overtopping  of  structures. 

Structure  crest  elevations  will  not  be  tested  or  optimized. 

No  wave-wave  or  wave-current  interaction. 

No  wave  breaking  effects. 

Diffraction  around  the  structure  ends  will  be  represented  by  diffraction 
around  a  blunt  vertical  wall  with  a  specified  reflection  coefficient. 

Within  limits  of  the  assumptions,  the  numerical  modeling  approach  can  be 
expected  to  give  a  reasonable  assessment  of  the  proposed  plans.  This  ap¬ 
proach  was  selected  because  POD’s  allowable  time  and  study  funds  were 
limited  and  design  alternatives  for  this  particular  project  were  extensive.  The 
procedure  of  this  study  is  described  in  the  following  paragraphs. 

The  deepwater  wave  conditions  for  the  southwest  coast  of  the  island  of 
Maui  were  established  from  the  Monitoring  Completed  Coastal  Projects 
(MCCP)  (US  Army  Corps  of  Engineers  (USACE)  1987)  measurements  taken 
at  Barbers  Point  Harbor,  Oahu,  HI.  The  percent  occurrences  of  the 
deepwater  conditions  were  calculated  to  later  determine  the  percent  occurrence 
of  the  wave  heights  inside  the  harbor.  The  method  to  establish  the  deepwater 
conditions  is  presented  in  Chapter  2  of  this  report,  "Deepwater  Wave 
Conditions." 

The  deepwater  waves  were  transformed  to  the  Maalaea  Harbor  vicinity 
through  application  of  the  SHALWV  numerical  model  (Hughes  and  Jensen 
1986).  SHALWV  simulates  growth,  decay,  propagation,  shoaling,  refraction, 
and  sheltering  of  a  directional  wave  spectrum  over  arbitrary  bathymetry.  The 
SHALWV  model  was  chosen  due  to  its  ability  to  include  effects  pertinent  to 
this  study  such  as:  (a)  simulation  over  an  extensive  area,  (b)  large  depth 
gradients  along  the  grid  boundaries,  (c)  input  along  two  grid  boundaries, 

(d)  wave  refraction  and  diffraction  around  islands,  and  (e)  the  ability  to  model 
refined  areas  through  subgrids.  The  SHALWV  model  is  presented  in 
Chapter  3  of  this  report,  "Wave  Transformation  Modeling." 

The  resulting  wave  conditions  of  SHALWV  were  then  used  as  input  to 
HARBD  to  determine  the  wave  response  inside  the  harbor.  HARBD  is  a 
steady  state  finite  element  model  that  calculates  linear  wave  oscillations  in 
harbors  of  arbitrary  configuration  and  variable  bathymetry.  The  effects  of 
bottom  friction  and  boundary  absorption  (reflection)  are  included.  Through 
application  of  HARBD,  the  resulting  wave  heights  in  the  harbor  entrance  and 
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access  channels,  turning  basin,  and  berthing  areas  were  determined  and  the 
percent  occurrence  of  those  conditions  were  calculated  using  the  results  of 
both  the  SHALWV  and  HARBD  models.  The  HARBD  model  and  details  and 
results  of  the  procedures  are  presented  in  Chapter  4  of  this  report,  "Wave 
Response  Modeling.” 
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2  Deepwater  Wave 
Conditions 


Maalaea  small  boat  harbor  is  affected  primarily  by  wave  conditions 
resulting  from  southern  swell  and  "Kona",  or  low  pressure  systems,  storm 
waves  arriving  from  south  or  south-southwest.  Sources  of  wave  information 
for  the  Hawaiian  Islands  and  particularly  for  use  at  the  Maalaea  Harbor 
vicinity  are  extremely  limited. 

Possible  sources  of  wave  information  in  the  Hawaiian  Islands  include 
(Figure  7);  the  Wave  Information  Studies  (WIS)  deepwater  hindcasts  for  the 
Pacific  coast  which  include  seven  stations  around  the  Hawaiian  Islands 
(Corson,  et  al.  1986);  the  deepwater  buoy.  Station  51(X)1  (Gilhousen,  et  al. 
1986),  operated  by  the  National  Oceanic  and  Atmospheric  Administration’s 
National  Data  Buoy  Center  (NDBC);  the  Summary  of  Synoptic  Meteorological 
Observations  (U.S.  Naval  Weather  Service  Command  1976)  climatological 
summaries  of  shipboard  wave  observations;  and  the  Monitoring  Completed 
Coastal  Projects  (MCCP)  Program  of  the  US  Army  Corps  of  Engineers 
(USACE  1987)  slope  array  at  Barbers  Point,  Oahu. 

Of  the  available  data  sources  in  the  Hawaiian  Islands,  the  wave 
measurements  from  the  MCCP  slope  array  at  Barbers  Point.  Oahu,  are  the 
only  directional  data  that  have  exposures  representative  of  the  coast  relevant  to 
Maalaea.  Although  Barbers  Point  is  on  a  different  island  than  Maalaea,  the 
coastline  orientation  and  southern  exposure  are  similar  at  the  two  sites. 


Deepwater  Wave  Climate 


Offshore  wave  climate  at  the  Maalaea  harbor  site  was  estimated  using  the 
MCCP  Barbers  Point  slope  array  measurements.  Although  data  sources  are 
limited,  the  percent  occurrence  of  significant  wave  height,  peak  spectral 
period,  and  direction  are  adequately  represented  in  the  Barbers  Point  data. 
These  data  include  effects  of  southern  swell,  which  is  an  important  factor 
relevant  to  the  Maalaea  harbor  site. 
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Figure  7.  Sources  of  wave  information  in  the  Hawaiian  Islands 


Maalaea  harbor  is  exposed  to  a  sector  of  wave  angles  from  about  the 
135.0-  to  270.0-deg  azimuth.  MCCP  measurements  within  this  sector  were 
taken  as  pan  of  the  deepwater  climate  at  Maalaea.  Measurements  outside  this 
sector  were  not  used  because  the  harbor  is  sheltered  by  land  from  those 
directions. 

The  MCCP  Barbers  Point  measurements  represent  a  general  wave  climate 
for  a  southwest-facing  coast.  Therefore,  measurements  taken  during  various 
intervals  over  the  period  from  July  1986  through  August  1990  were  compiled 
to  obtain  a  data  set  representative  of  one  complete  year.  The  measurements 
were  then  "unrefracted"  from  28  ft  to  an  approximate  1,300-ft  depth  to 
transfer  the  data  set  to  an  appropriate  location  fronting  Maalaea.  This  process 
was  accomplished  through  application  of  an  Automated  Coastal  Engineering 
System  (ACES)  code  which  implements  Snell’s  Law.  Bottom  contours  were 
assumed  straight  and  parallel.  The  percent  occurrences  from  each  direction  in 
the  exposed  sector  were  then  compiled  and  are  given  in  Appendix  A  and 
summarized  by  direction  in  Table  1  at  the  end  of  the  text  of  this  report.  The 
compiled  percent  occurrence  of  height,  period,  and  direction  are  considered  as 
the  best  possible  representation  of  the  deepwater  wave  climate  for  Maalaea 
Harbor.  The  percent  occurrence  of  waves  corresponding  to  Table  1  is  shown 
in  Figure  8. 
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Figure  8.  Percent  occurrence  of  deepwater  wave  conditions  by  direction 
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3  Wave  Transformation 
Modeling 


Wave  Transformation  Model 


Once  the  deepwater  wave  conditions  are  established,  it  is  necessary  to 
transform  the  waves  shoreward  to  the  Maaiaea  vicinity.  This  was  accom¬ 
plished  using  the  numerical  wave  model  SHALWV  (Jensen.  Vincent,  and 
Abel  1987).  SHALWV  numerically  simulates  growth,  decay,  propagation, 
shoaling,  reft-action,  and  sheltering  of  a  directional  wave  spectrum  over 
arbitrary  bathymetry.  The  modeled  spectra  are  represented  as  fully  two- 
dimensional  spectra  in  discrete  frequency  and  direction  bands.  The  model  is 
time-dependent  and  is  most  often  used  to  simulate  time-varying  wave  and 
wind  conditions.  SHALWV  was  chosen  to  transform  waves  to  the  Maaiaea 
vicinity  due  to  the  sizable  extent  of  the  grid  required  to  include  the  island  of 
Kahoolawe.  For  this  application,  winds  are  not  used:  therefore,  no  wave 
growth  occurs  and  wave  input  is  held  constant  for  each  wave  condition.  The 
fetch  of  the  modeled  area  is  approximately  20  miles,  therefore  the  effect  of 
local  winds  is  limited  to  higher  frequencies.  Although  some  differences  in  the 
distribution  of  wave  energy  in  frequency  and  direction  can  be  expected  with 
the  introduction  of  local  winds,  only  extreme  conditions  will  contribute 
significantly.  Holding  the  wave  conditions  constant,  in  this  case  for  5  hr, 
neglects  changes  in  regional  winds  and  waves  on  a  smaller  time  scale.  This  is 
a  limitation  of  the  measured  data,  which  was  collected  in  6-hr  increments. 

The  model  is  based  on  the  solution  of  the  inhomogeneous  energy  balance 
equation  solved  with  finite-difference  methods  using  square  grid  cells  to 
describe  the  bathymetry  and  wind  field.  The  field  equation  represents  wind 
and  wave  growth,  refraction,  shoaling,  nonlinear  wave-wave  interactions,  high 
frequency  energy  dissipation,  wave  bottom  interactions,  and  decomposition  of 
the  energy  into  wind-sea  and  swell  wave  components.  The  maximum  time- 
step  interval  is  determined  from  the  following  Courant  number  stability 
criterion  which  is  based  on  grid  size  and  water  depth; 

—  >  CJf) 

A/  ^ 
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where 


tL  =  length  of  grid  cell 
tj  =  computational  time  step 

Cg  =  group  velocity  associated  with  lowest  frequency  at  the 
deepest  grid  point 
f  —  lowest  spectral  frequency 

The  Courant  number  criterion  insures  that  wave  energy  does  not  propagate 
more  than  one  grid  cell  during  a  time  step  (Hughes  and  Jensen  1986). 


Wave  Transformation  Simulation 


The  SHALWV  model  uses  a  rectangular,  uniformly  spaced,  finite  differ¬ 
ence  grid  that  extends  over  the  general  area  of  interest.  If  necessary,  the 
model  allows  additional  subgrids  that  are  used  to  resolve  complex  bathymetry 
in  areas  of  interest.  The  main  grid  used  in  this  study,  shown  in  Figure  9,  has 
50  cells  alongshore  (positive  x-axis  directed  east)  and  42  cells  across-shore 
(positive  y-axis  directed  north).  The  grid  has  square  cells  of  about  3,000  ft  in 
the  X  and  y  directions.  This  grid  covered  the  region  of  interest  including 
Maalaea  Bay  and  important  sheltering  features  such  as  the  island  of 
Kahoolawe  and  McGregor  Point.  Waves  were  input  along  the  x  or  y  axis 
depending  on  the  direction  being  modeled. 

In  order  to  resolve  the  bathymetry  near  Maalaea  Harbor,  a  higher  resolu¬ 
tion  subgrid  of  Maalaea  Bay  was  included  in  the  modeling  effort.  The  extent 
of  the  subgrid  is  outlined  in  Figure  9.  The  subgrid  consists  of  35  alongshore 
and  31  across-shore  square  grid  cells  with  about  750-ft  sides.  The  resolution 
of  the  subgrid  is  approximately  four  times  greater  than  the  main  grid.  Wave 
output  at  the  interface  between  the  subgrid  and  the  HARBD  outer  boundary 
were  used  as  input  to  HARBD. 

SHALWV  transformation  estimates  were  performed  from  a  depth  of 
approximately  1,300  ft  at  the  offshore  boundary  to  approximately  30  ft  at  the 
HARBD  outer  boundary  (Maalaea  harbor).  Representative  period-direction 
combinations  were  selected  from  the  modified  slope  array  data  for  input  to 
SHALWV.  Direction  bands  modeled  were  22.5  deg  and  centered  about  the 
135.0-,  157. 5-,  180.0-,  202. 5-,  225.0-,  247.5-,  and  270.0-deg  azimuth,  as 
shown  in  Figure  9.  Peak  periods  ranged  from  9  to  20  sec  and  selected  wave 
heights  ranged  from  3  to  8  ft  based  on  shoaling  estimates  which  would  exceed 
the  maximum  I-  and  2-ft  wave  height  criterion  discussed  in  Chapter  4. 

The  wave  conditions  input  to  SHALWV  and  the  corresponding  transformed 
wave  height  and  direction  are  given  in  Table  2  at  the  end  of  this  text.  Figures 
10  through  16  are  plots  of  wave  direction  vectors  that  show  representative 
wave  refraction  from  the  SHALWV  main  grid  input  boundary,  (x  or  y  axis), 
to  the  subgrid  boundary  in  the  Maalaea  Harbor  vicinity.  Each  figure  shows 
resulting  wave  direction  vectors  from  input  conditions  of  an  8-ft  wave  height. 
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Figure  9.  Maalaea  SHALWV  grid  and  bathymetry 
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Figure  10.  Wave  direction  vectors  for  135.0  deg 
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Figure  1 1 .  Wave  direction  vectors  for  1 57.5  deg 


Inpul :  Dir=202 1i=20s  H=8  fl  Inpul :  Dir=202  T=lls  H=8  fl 
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direction  vectors 


Figure  14.  Wave  direction  vectors  for  225.0  deg 


direction  vectors 


Inpul :  Dir=e70  'P=20s  H=8  ft  Input :  Dir=270  T=lls  H=8  ft 

Wave  DirecUon  Vectors  Wave  Direction  Vectors 


Figure  16.  Wave  direction  vectors  for  270.0  deg 


an  11-  and  20-sec  wave  period,  and  wave  directions  ranging  from  the  135.0- 
to  270.0-deg  azimuth.  The  length  and  direction  of  the  vectors  represent  the 
relative  wave  height  and  direction,  respectively,  for  each  grid  point. 

These  plots  show  that  waves  coming  from  247.5,  225.0.  202.5.  and 
180.0  deg  produce  the  most  wave  energy  at  the  Maalaea  Harbor  vicinity. 

Also,  the  shorter  wave  periods  (11  sec)  result  in  a  higher  percentage  of  energy 
reaching  the  harbor  area  than  the  longer  period  (20-sec)  waves.  Waves 
coming  from  247.5  and  225.0  deg  travel  in  nearly  straight  line  paths  to  the 
harbor  area,  resulting  in  a  shoaling  dominant  transformation  process.  Waves 
from  202.5  deg  begin  refracting  toward  the  island  of  Kahoolawe  on  the  east 
and  west  sides,  which  results  in  wave  refraction  focusing  toward  the  harbor 
vicinity.  Note  that  the  eftects  of  wave  diffraction  around  the  island  of 
Kahoolawe  are  not  included  in  the  simulation.  These  effects  do  not  signifi¬ 
cantly  affect  the  wave  energy  reaching  the  harbor  area  for  wave  periods  of 
20  sec  or  less,  since  the  distance  between  the  island  and  the  harbor  area 
(12.5  miles)  is  much  greater  than  20  wavelengths  (7.8  miles  for  a  20-sec 
wave).  Waves  coming  from  180.0  deg  travel  through  Alalakeiki  Channel  and 
are  parallel  to  the  bottom  contours  in  Maalaea  Bay.  This  results  in  another 
direct  approach  of  wave  energy  into  the  harbor  vicinity. 

Given  the  percent  occurrence  of  waves  from  the  directions  between  135.0 
and  270.0  deg,  (Figure  8),  the  most  significant  wave  directions  are  247.5 
(16.37  percent),  225.0  (50.78  percent),  and  202.5  (19.9  percent)  deg.  The 
percent  occurrence  of  waves  coming  from  180.0  deg  is  2.2  percent,  and  from 
157.5  and  135.0  deg  is  less  than  I  percent.  Therefore,  the  frequency  of 
waves  affecting  the  harbor  from  these  directions  will  be  minimal. 
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4  Harbor  Wave  Response 
Modeling 


Harbor  Wave  Response  Model 


The  numerical  model  HARBD  (Chen  and  Houston  1987)  was  used  to 
simulate  the  wave  response  at  Maalaea  small  boat  harbor,  Maui,  Hawaii. 
HARBD  is  a  steady  state  hybrid  finite  element  model  that  calculates  linear 
wave  oscillations  in  harbors  of  arbitrary  configuration  and  variable 
bathymetry.  The  model  is  advantageous  over  other  numerical  harbor  models 
since  bottom  friction  and  boundary  reflection  are  included.  The  bottom  fric¬ 
tion  is  assumed  to  be  proportional  to  flow  velocity  with  a  phase  difference. 
The  boundary  absorption  is  based  on  a  formulation  similar  to  that  in  the 
impedance  condition  in  acoustics  and  is  expressed  in  terms  of  wave  number 
(2ir/L  where  L  is  the  wavelength)  and  reflection  coefficient  of  the  boundary. 
The  result  is  that  HARBD  predicts  wave  amplitudes  that  are  more  realistic 
than  those  from  previous  models  (Chen  and  Houston  1987).  HARBD  was 
originally  developed  for  harbor  oscillations  (long-period  waves),  and  the 
general  formulation  was  adapted  for  wind  waves  (short-period  waves)  by 
Houston  (1981). 

The  model  has  been  tested  and  compared  with  known  analytical  solutions 
for  a  number  of  cases  and  the  results  are  excellent  (Chen  1984,  and  Chen  and 
Houston  1987).  It  has  been  applied  in  the  design  or  modification  of  Agat 
Harbor,  Guam  (Farrar  and  Chen  1987);  Kawaihae  Harbor.  Hawaii.  Hawaii 
(Lillycrop,  Bratos,  and  Thompson  1990):  and  Barbers  Point  Harbor,  Oahu. 
Hawaii  (Durham  1978).  The  model  was  instrumental  in  studying  the  effects 
of  entrance  channel  dredging  at  Morro  Bay  Harbor,  California  (Kaihatu, 
Lillycrop,  and  Thompson  1989),  and  analyzing  harbor  resonance  at  Los 
Angeles-Long  Beach  Harbor  California  (Sargent  1989).  The  model  was  used 
to  plan  wave  protection  at  Fisherman’s  Wharf.  San  Francisco,  California 
(Bottin,  Sargent,  and  Mize  1985);  Green  Harbor,  Massachusetts  (Weishar  and 
Aubrey  1986);  Los  Angeles-Long  Beach  Harbor.  California  (Houston  1976); 
and  to  estimate  the  wave  conditions  in  Indiana  Harbor,  Indiana  during  a  study 
of  sediment  disposal  alternatives  (Clausner  and  Abel  1986).  HARBD  was 
compared  to  laboratory  data  collected  from  the  physical  model  study  of 
Barcelona  Harbor,  Buffalo,  New  York  (Crawford  and  Chen  1988)  with 
encouraging  results.  The  predictions  of  HARBD  are  currently  under  further 
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comparison  with  prototype  and  physical  model  data  collected  from  recent 
studies  of  Barbers  Point  Harbor.  Oahu.  Hawaii. 

HARBD  uses  a  hybrid  element  method  in  which  a  finite  element  solution 
in  the  interior  region  of  the  harbor  is  matched  to  an  analytical  solution  in  the 
exterior  region.  In  the  interior  region,  HARBD  allows  arbitrary  bathymetry, 
(i.e.,  shallow,  intermediate,  and  deepwater  waves),  variable  configuration, 
and  the  effects  of  bonom  friction  and  boundary  reflection. 

In  model  formulation  for  arbitrary  depth  water  waves,  the  water  domain  is 
divided  into  near  and  semi-infinite  far  regions.  The  near  region  is  bounded  by 
an  artificial  180  deg  semi-circular  boundary  outside  the  harbor  and  includes 
the  interior  harbor  and  all  marine  structures  and  bathymetry  of  interest.  The 
far  region  is  an  infinite  semicircular  ring  shape  bounded  by  the  near  region 
and  the  coastlines.  The  region  extends  to  infinity  in  all  horizontal  directions. 
The  semi-infinite  far  region  is  assumed  to  have  a  constant  water  depth  and  no 
bottom  friction  (Chen  and  Houston  1987).  The  finite  near  region,  which 
contains  the  area  of  interest,  is  subdivided  into  a  mesh  of  triangular  shaped 
elements.  The  length  of  the  sides  of  each  element  is  determined  from  the 
desired  grid  resolution  and  design  wave  parameters.  The  water  depth  and 
bottom  friction  coefficient  are  specified  for  each  element,  and  a  reflection 
coefficient  is  assigned  to  each  element  on  the  solid  boundaries.  The  model 
requires  a  wave  period  and  direction  as  input.  The  solution  consists  of  an 
amplification  factor  (i.e.,  the  ratio  of  the  wave  height  to  the  incident  wave 
height)  and  a  corresponding  phase  angle  for  the  entire  near  region.  The  phase 
angle  is  of  little  imponance  to  the  present  study. 

The  governing  partial  differential  equation  is  derived  through  application  of 
linear  wave  theory  to  the  continuity  and  momentum  equations.  This  also 
assumes  all  dependent  variables  are  periodic  in  time  with  angular  frequency  u. 
These  steps  yield  the  following  generalized  Helmholtz  equation  (Chen  1986) 
in  which  the  velocity  potential  <(>  is  solved: 

V<Xcc  V<t>)  +  =  0  (2) 

'  c 

where 

V  =  horizontal  gradient  operator 

X  =  complex  bottom  friction  factor 

c  =  wave  phase  velocity  =  («/*) 

Cg  =  wave  group  velocity  =  [c/2{l  (2xh/sinh  2ich)}l 

K  =  wave  number,  (2ir/L),  where  L  =  wavelength 
w  =  angular  frequency 

<f>  =  velocity  potential 

h  =  water  depth 
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The  wave  number  is  obtained  from  the  dispersion  relation. 


^  tanh  (kA) 

where  g  =  acceleration  due  to  gravity 

The  complex  bottom  friction  factor  X  is  assumed  proportional  to  the 
maximum  velocity  at  the  bottom  and  is  defined  as: 


X 


1 


1  + 


h  sinh  ich 


expfiy) 


(4) 


where 

/3  =  dimensionless  bottom  friction  coefficient  that  varies  spatially 

a^  =  incident  wave  amplitude 

7  =  phase  shift  between  stress  and  flow  velocity 

i  =  (-1)*'= 


The  effects  of  bottom  friction  do  not  necessarily  need  to  be  included  in  the 
general  solution.  This  is  accomplished  by  setting  /3  =  0,  which  results  in 
X  =  1,  and  Equation  1  reduces  to  an  expression  which  excludes  bottom 
friction. 


For  the  absorptive  boundary  condition  along  the  solid  harbor  boundaries, 
the  model  adopts  the  impedance  condition  used  in  acoustics  in  terms  of  the 
boundary  reflection  coefficient  expressed  as: 


-  a4>  =  0 
dn 

with 


a 


IK 


1  -K 
1  . 


(5) 


(6) 


where 

a  =  dimensional  coefficient  related  to  the  boundary  reflection 
n  =  unit-normal  vector  directed  outward  from  the  fluid  domain 


Similar  to  the  friction  coefficient,  when  K,  =  1 ,  then  a  =  0  and  Equation  5 
reduces  to  a  zero  velocity  potential  normal  to  the  boundary  (Sargent  1989). 
This  infers  a  perfectly  reflecting  boundary  condition. 

A  conventional  finite  element  approximation  is  used  in  the  near  region,  and 
an  analytical  solution  with  unknown  coefficients  is  used  to  describe  the  semi- 
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infinite  far  region.  Conditions  in  the  near  and  far  regions  must  be  matched 
along  the  artificial  semicircle  boundary.  This  requirement  is  met  by  HARBD 
routines  which  automatically  match  the  solutions  using  the  stationarity  of  a 
functional,  to  a  series  of  Hankel  Functions  which  give  the  solution  for  the 
infinite  region  (Farrar  and  Chen  1987).  The  hybrid  element  numerical  tech¬ 
niques  used  in  the  formulation  are  discussed  in  greater  detail  in  Chen  and  Mei 
(1974). 

The  HARBD  model  is  intended  to  simulate  waves  that  can  be  adequately 
described  by  the  mild  slope  equation  (Equation  2).  Model  accuracy  decreases 
as  wave  conditions  approach  those  outside  the  validity  of  this  governing 
equation.  HARBD  does  not  simulate  nonlinear  processes  such  as  wave 
breaking,  wave  transformation  and  overtopping  of  structures,  and  wave 
current  interaction;  however,  the  model  pr^icts  wave  heights  accurately  if 
these  processes  are  not  dominant.  Since  nonlinear  processes  naturally  occur 
in  the  prototype,  care  and  consideration  of  the  effects  must  be  taken  in 
interpretation  of  results. 


Finite  Element  Grids 


Finite  element  grids  generated  for  the  Existing  Plan  and  Plans  1,2,  and  3 
are  shown  in  Figures  17  through  20,  respectively.  The  grid  used  for  both 
Plans  la  and  lb  is  shown  in  Figure  21.  All  grids  cover  approximately  the 
same  interior  harbor  areas,  however,  coverage  of  the  offshore  area  is  variable. 
In  order  to  model  areas  pertinent  to  each  plan,  the  Existing  and  Plans  1.  la. 
and  lb  include  the  offshore  area  eastward,  while  Plans  2  and  3  include  the 
offshore  area  westward.  The  radius  of  the  semicircular  boundary  is  approxi¬ 
mately  800  ft  for  the  existing  and  Plan  1.1a.  and  lb  grids,  and  approximately 
700  ft  for  the  Plan  2  and  3  grids.  The  radial  distance  is  designed  to  include 
the  entrance  channel  and  allow  enough  area  to  include  possible  modifications. 

Total  numbers  of  elements  (triangles),  nodes  (triangular  corners),  and 
boundary  elements  are: 

Existing;  7,146  elements.  3,752  nodes.  252  boundary  elements 

Plan  1:  6,765  elements,  3,613  nodes.  356  boundary  elements 

Plan  2:  7,866  elements,  4,176  nodes,  353  boundary  elements 

Plan  3:  7,911  elements,  4,215  nodes,  386  boundary  elements 

Plans  la,  lb:  6,810  elements,  3,636  nodes,  357  boundary  elements 

Each  grid  was  designed  with  a  grid  resolution  of  approximately  six  elements 
per  wavelength,  based  on  an  8-sec  wave  period  and  a  basin  depth  of  8  ft. 
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Figure  17.  Finite  element  grid  for  Existing  Plan 


Figure  1 8.  Finite  element  grid  for  Plan  1 
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Figure  19.  Finite  element  grid  for  Plan  2 


Figure  20.  Finite  element  grid  for  Plan  3 
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Figure  21 .  Finite  element  grid  for  Plans  la  and  1b 

The  grid  bathymetry  was  obtained  from  POD  hydrographic  surveys  taken 
in  1989.  The  design  still-water  level  (swl)  was  us^  for  all  tests.  Boundary 
reflection  coefficients  of  the  interior  harbor  walls  were  calculated  using  meth¬ 
ods  in  the  Shore  Protection  Manual  (SPM  1984).  and  were  refined  upon  rec¬ 
ommendation  from  POD.  The  assigned  reflection  coefficients  for  the  existing 
structures  were  0.40  for  the  1 :2-sloped  seaward  and  harbor  sides  of  the  south 
breakwater  and  along  the  west  and  north  walls  of  the  harbor  basin.  1 .0  for  the 
paved  wharf,  and  0.35  for  the  l;2-sloped  seaward  and  harbor  sides  of  the  east 
breakwaters.  For  Plan  1,  the  reflection  coefficients  for  the  structural  improve¬ 
ments  were  0.25  for  the  1  ;2-sloped  south  breakwater  extension  and  seaward 
side  revetment  and  0.35  for  the  1 :!  .5-sloped  interior  revetments  along  the 
existing  east  breakwater  and  at  the  center  of  the  harbor  basin.  For  the  struc¬ 
tural  improvements  of  Plan  2,  the  reflection  coefficients  were  0.25  for  the 
l:2-sloped  east  breakwater  extension  and  0.35  for  the  1:1. 5-sloped  interior 
revetments  as  in  Plan  1.  The  additional  250-ft-long  l:2-sloped  east  break¬ 
water  extension  in  Plan  3  was  assigned  a  reflection  coefficient  of  0.25. 
Reflection  coefficients  in  Plan  la  were  the  same  as  in  Plan  1.  Plan  lb 
differed  from  Plan  la  in  that  the  reflection  coefficient  along  the  east  face  of 
the  center  mole  was  set  to  1 .00.  The  open  boundary  along  the  east  diameter 
of  the  semicircle  was  fully  transmissive.  The  bottom  friction  factor  /3  was  set 
at  0.05  since  the  entire  bottom  was  sandy  (Kaihatu,  Lillycrop,  and  Thompson 
1989). 

The  grids  for  this  study  were  generated  through  application  of  automated 
finite  element  grid  generation  software  developed  at  the  Oregon  Graduate 
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Institute  (OGI)  by  Dr.  Antonio  M.  Baptista  and  Mr.  Paul  J.  Turner.  This 
procedure  for  automated  HARBD  grid  generation  was  an  original  WES 
application. 

Harbor  Wave  Response  Simulation 


To  establish  the  wave  climate  incident  to  Maalaea  harbor,  a  total  of  187 
deepwater  wave  height,  period,  and  direction  combinations  were  input  to  the 
SHALWV  model  and  transformed  to  the  Maalaea  harbor  vicinity.  The 
selected  SHALWV  output  locations  and  their  relation  to  Maalaea  Harbor  were 
given  previously  in  Chapter  3  of  this  report.  Deepwater  wave  characteristics 
and  the  resulting  transformed  wave  conditions  are  given  in  Table  2.  To 
determine  wave  heights  throughout  the  harbor,  the  resulting  SHALWV  wave 
heights  were  multiplied  with  the  HARBD  amplification  factors  corresponding 
to  each  deepwater  condition.  The  187  wave  height,  period,  and  direction 
combinations  were  tested  for  the  Existing  Plan  and  Plans  1.  2.  3,  la,  and  Ib. 
All  simulations  were  run  on  the  WES  CRAY  Y-MP  supercomputing  facilities. 

Output  "basins"  were  selected  for  each  plan  tested  to  determine  wave 
response  throughout  the  harbor.  A  basin  is  an  area  consisting  of  a  specified 
number  of  elements  from  which  the  mean  value  of  the  results  of  those 
elements  is  calculated.  Sixteen  output  basin  locations  were  selected  for  the 
Existing  Plan,  23  for  Plans  1,  2,  ia,  and  lb,  and  24  for  Plan  3.  The  locations 
were  selected  in  areas  of  interest  for  safe  navigation  and  mooring  by  CERC 
and  POD,  and  are  shown  for  each  plan  tested  in  Figures  22  through  25.  For 
the  Existing  and  Plans  1,  2,  and  3,  basins  1  through  3,  1  through  5,  1  through 
6,  and  2  through  7,  respectively,  are  located  throughout  the  harbor  entrance 
and  access  channels  and  turning  basin  with  a  2-ft  maximum  wave  height 
criterion.  Basins  3  through  16,  6  through  23,  7  through  23.  and  8  through 
24,  respectively,  are  located  in  the  harbor  berthing  areas  with  a  1  ft  maximum 
wave  height  criterion.  Basin  locations  and  numbers  in  Plans  la  and  Ib  are 
identical  to  those  in  Plan  1.  The  resulting  HARBD  amplification  factors  at 
these  basiho  lor  each  deepwater  wave  condition  were  saved  and  tabulated  for 
each  plan  (Tables  3  through  44). 

The  percent  occurrence  of  wave  heights  exceeding  1  ft  in  the  berthing 
areas  and  2  ft  in  the  entrance  and  access  channels  and  turning  basin  were 
calculated  for  the  Existing  and  Plans  1,  2,  3.  la,  and  lb.  The  procedure  to 
calculate  the  percent  occurrence  of  wave  heights  exceeding  the  I -ft  maximum 
criterion  is  as  follows.  The  largest  HARBD  amplification  factor  of  the  basins 
located  in  the  1-ft  maximum  wave  height  criterion  areas  (berthing  areas),  was 
selected  for  each  deepwater  wave  condition.  The  selected  HARBD  amplifica¬ 
tion  factors  were  then  multiplied  by  the  transformed  wave  heights  from 
SHALWV  corresponding  to  each  deepwater  wave  period  and  direction  with  a 
wave  height  of  3  ft.  If  the  resulting  SHALWV-HARBD  wave  height  does  not 
exceed  1  ft,  the  largest  HARBD  amplification  was  then  multiplied  by  the 
transformed  SHALWV  wave  height  corresponding  to  the  deepwater  wave 
conditions  with  a  wave  height  of  4  ft.  The  iterative  process  contin.  ed  until 
the  resulting  wave  height  exceeded  the  maximum  1  ft  wave  height  criterion  or 
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Figure  24.  Output  basin  locations  for  Plan  2 


-igure  25.  Output  basin  locations  for  Plan  3 
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the  deepwater  wave  height  exceeded  9  ft.  The  percent  occurrence  of  wave 
heights  exceeding  9  ft  was  included  in  all  tabulations  since  the  incremental 
deepwater  wave  heights  greater  than  9  ft  could  not  be  defined.  The  percent 
occurrence  of  those  resulting  wave  heights  which  exceeded  I  ft  were  then 
tabulated  from  the  percent  occurrence  tables  for  the  deepwater  wave  condi¬ 
tions.  The  same  procedure  was  used  to  calculate  the  percent  occurrence  of 
wave  heights  exceeding  the  maximum  2-ft  criterion  for  those  basins  located  in 
the  harbor  entrance  and  access  channels  and  turning  basin.  As  expected,  wave 
heights  from  the  Existing  Plan  are  largest.  Wave  heights  from  Plan  2  exceed 
those  of  Plans  1  and  3,  and  wave  heights  from  Plan  3  are  lowest  overall. 

The  above-mentioned  procedure  can  be  followed  to  calculate  the  resulting 
HARBD-SHALWV  wave  height  at  a  specific  output  location  of  the  Existing  or 
Plans  1,  2,  3,  la,  and  lb  for  a  specified  deepwater  wave  period,  direction, 
and  height.  Tabulations  of  these  results  are  not  included  in  the  text  due  to  the 
substantial  amount  of  data  involved. 

Tables  45  through  50  are  tabulations  of  the  HARBD-SHALWV  wave 
heights  initially  exceeding  the  HQUSACE  criterion  for  each  deepwater  wave 
direction.  The  deepwater  wave  period,  exceeding  wave  height,  deepwater 
wave  height,  HARBD  amplification  factor,  SHALWV  wave  height,  and  basin 
in  which  they  occurred  are  given  for  the  1-  and  2-ft  maximum  criteria  for  the 
Existing  and  Plans  1,  2,  3,  la,  and  lb.  For  the  Existing  Plan,  Table  45 
shows  that  the  wave  heights  initially  exceeding  the  maximum  1  ft  criterion 
berthing  areas  (basins  3  through  16)  were  caused  by  a  9-sec  wave  from 

247.5  deg,  a  13-sec  wave  from  directions  of  157.5  and  225.0  deg,  a  15-sec 
wave  from  202.5  deg,  and  a  17-sec  wave  from  180.0  deg.  The  exceeding 
wave  heights  from  directions  of  157.5,  180.0,  202.5,  and  247.5  deg  occurred 
in  basin  4  and  the  wave  from  225.0  deg  occurred  in  basin  9.  Referring  to 
Figure  22,  basin  4  is  located  just  inside  the  entrance  toward  the  eastern  side  of 
the  harbor  and  basin  9  is  located  along  the  northern  wall.  Wave  heights 
exceeding  the  2-ft  maximum  criterion  turning  basin  and  entrance  channel 
(basins  1  and  2)  resulted  from  a  9-sec  wave  from  directions  of  225.0  and 

247.5  deg,  a  15-sec  wave  from  180.0  deg,  and  a  17-sec  wave  from 

202.5  deg.  These  waves  occurred  at  the  harbor  entrance  in  basin  I. 

In  evaluating  Table  46  for  Plan  1,  the  wave  heights  initially  exceeding  the 
l-ft  criterion  in  the  berthing  areas  (basins  6  through  23)  are  caused  by  a 
13-sec  wave  from  180.0  deg  and  a  15-sec  wave  from  225.0  deg.  The  wave 
from  180.0  deg  occurred  in  basin  11,  located  along  the  existing  wharf,  and 
the  225.0-deg  wave  occurred  in  basin  8.  located  east  of  the  interior  revetment. 
The  maximum  2-ft  wave  height  criterion  was  not  exceeded  for  Plan  1 . 

Wave  conditions  initially  exceeding  the  maximum  l-ft  criterion  berthing 
areas  for  Plan  2  (basins  7  through  23)  include:  an  1 1  -sec  wave  from 

157.5  deg,  occurring  along  the  south  breakwater  in  basin  23;  a  20-sec  wave 
from  180.0  deg  in  basin  8;  and  a  20-sec  wave  from  202.5  deg  and  an  1 1-sec 
wave  from  225.0  deg,  both  occurring  in  basin  7.  Basins  7  and  8  are  located 
near  the  east  breakwater.  Wave  heights  initially  exceeding  the  2-ft  maximum 
criterion  channels  and  turning  basin  include:  an  1 1-sec  wave  from  157.5  deg; 


36 


ChaptAr  4  Harbor  Wave  Raaponae  Modeling 


a  20-sec  wave  from  180.0  deg  and  202.5  deg;  and  a  9-sec  wave  from 
225.0  deg.  These  waves  all  occurred  in  basin  1  located  at  the  harbor 
entrance. 

As  shown  in  Table  48,  none  of  the  deepwater  wave  conditions  resulted  in 
wave  heights  exceeding  the  maximum  1-  and  2-ft  criteria  for  Plan  3,  however, 
the  percent  occurrence  of  wave  heights  greater  than  9  ft  was  included  in  the 
tabulations  for  this  plan. 

For  Plan  la,  wave  heights  exceeded  the  I -ft  berthing  area  criterion  for 
many  of  the  shorter  wave  periods  with  directions  from  180  deg  to  247.5  deg 
(Table  49).  In  most  cases,  the  exceedance  occurred  at  basin  1 1,  located  along 
the  existing  wharf.  The  9-sec  waves  caused  exceedances  at  basins  6  and  8, 
located  near  the  east  breakwater.  The  2-tt  criterion  was  exceeded  at  basin  1 
for  9  sec  periods  from  202.5  deg  and  225.0  deg  and  for  1 1-sec  periods  from 
225.0  deg. 

The  1-ft  criterion  was  exceeded  for  Plan  lb  for  many  of  the  shorter  periods 
with  directions  between  157.5  deg  and  247.5  deg  (Table  50).  The 
exceedances  generally  occurred  at  basin  11,  though  basins  6  and  7  also 
appeared.  The  2-ft  criterion  was  exceeded  at  basin  1  for  9-sec  waves  from 
202.5  deg  and  225.0  deg  and  1 1-sec  waves  from  225.0  deg. 

The  percent  occurrence  of  wave  heights  exceeding  the  maximum  1-  and 
2-ft  criteria  more  than  approximately  10  percent  of  the  time  were  calculated 
using  the  percent  occurrence  tables  of  deepwater  conditions  and  HARBD- 
SHALWV  wave  height  results  for  all  plans.  These  results  are  given  in  Tables 
51  through  62  and  illustrated  in  Figures  26  through  31.  Although  wave 
breaking  was  not  taken  into  account  in  Tables  51  through  62,  the  higher  wave 
heights  would  most  likely  have  broken  over  the  reef,  thus  reducing  wave 
heights  in  the  harbor.  In  evaluating  the  resulting  percent  occurrence  tables 
(Tables  51  through  62)  and  Figures  26  through  31,  it  is  apparent  that  waves 
approaching  from  the  southeast  (135.0-  to  157.0-deg)  directions  are  insignifi¬ 
cant  in  comparison  to  waves  approaching  from  the  south  to  west  (180.0-  to 
270.0-deg)  directions. 

The  percentage  of  wave  heights  exceeding  the  maximum  I  ft  and  2  ft 
criteria  for  the  Existing  and  Plans  1,  2,  3,  la.  and  lb  are  summarized  in 
Table  63,  along  with  the  HQUSACE  criteria.  These  values  are  conservative 
since  they  represent  basins  with  the  largest  wave  heights  occurring  in  the 
harbor  for  each  deepwater  wave  condition.  The  Existing  Plan  and  Plans  2  and 
lb  allow  one  or  both  of  the  HQUSACE  criteria  of  wave  heights  greater  than  1 
and  2  ft  more  than  10  percent  of  the  time  per  year  to  be  exceeded.  However, 
Plans  1  and  la,  which  include  structural  modification  to  the  east,  and  Plan  3, 
which  includes  structural  modification  toward  the  west,  satisfy  the  HQUSACE 
criteria  for  providing  adequate  protection  inside  the  harbor.  Although  Plan  la 
satisfies  the  1  fr  berthing  area  criterion,  the  results  indicate  that  basin  1 1 . 
along  the  existing  wharf,  may  be  marginally  acceptable. 
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5  Conclusions 


The  numerical  modd  studies  and  results  described  in  diis  report  should  be 
semi  in  light  of  die  following  considerations: 

a.  The  deqiwater  waves  were  based  on  measurements  from  MCCP  data 
collected  at  Barbers  Point,  Oahu.  Availability  of  incident  wave  data  at 
the  Maalaea  Harbor  vicinity  would  improve  foe  validity  of  foe  overall 
results. 

b.  The  revised  SPM  reflection  coefficients  were  based  on  estimates  from 
POD  and  were  not  reevaluated.  Researdi  in  this  area  continues  at 
CERC  for  bettCT  guidance. 

c.  The  following  assunqitions  were  made  in  foe  in^ilementation  of  foe 
HARBD  numerical  r^d  used  in  fois  study:  (1)  The  modd  does  not 
consider  wave  transmission  dirough  foe  breakwater,  overtopping  of 
structures,  and  wave  breaking  effects  in  foe  entrance  diannel. 

(2)  Structure  crest  devations  were  not  tested  or  optimized. 

(3)  Currents  and  nonlinear  effects  were  neglected.  (4)  Diffraction 
around  foe  structure  ends  was  represented  by  diffradion  around  a 
blunt  vertical  wall  wifo  specified  reflection  coefficients.  If  wave 
transmission  through  foe  breakwater  and  overtopping  of  structures  did 
occur  in  foe  harbor,  foe  increased  energy  would  result  in  iargd  wave 
heights  than  predicted.  The  presence  of  wave  currmits  and  breaking 
would  increase  hazardous  navigation;  howevd,  wave  breaking  would 
reduce  foe  energy  in  foe  harbor  and  result  in  lower  wave  heists  than 
predicted.  The  primary  effects  which  must  be  considered  within  a 
harbor  such  as  Maalaea  are  wave  refraction,  diffraction,  and 
dissipation  effects,  for  which  foe  model  has  been  well  vmifled. 

d.  The  HARBD  model  uses  monochromatic  waves  only. 

e.  The  resulting  percent  occurrence  of  wave  heights  exceeding  foe  1-  and 
2-ft  criteria  is  based  on  incidmtt  significant  wave  heists  H1/3.  The 
use  of  H|o  or  Hj  would  increase  foe  incident  wave  heights  by 
iq>proxiniately  27  and  67  percent,  respectively.  Therefore,  resulting 
wave  heists  inside  foe  harbor  would  increase,  and  foe  percent 
occurrence  of  wave  heists  exceeding  foe  1-  and  2-ft  wave  hei^t 
criteria  would  increase. 
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Based  on  the  results  of  this  study,  the  following  conclusions  were  reached: 

a.  The  Existing  Plan  is  unsatisfactory  in  providbg  the  harbor  with  ade¬ 
quate  protection  against  the  incident  wave  climate. 

b.  The  POD  plan  based  on  the  GDM  (Plan  1),  which  is  directed  toward 
the  east,  is  satisfactory  relative  to  &e  HQUSACE  design  criteria  for 
protecting  the  harbor  from  the  incident  wave  climate.  Plan  1  is  a 
recommmided  alternative. 

c.  Plan  2,  whidi  is  directed  toward  the  west,  will  not  protect  the  harbor 
adequatdy  from  Ae  deqpwater  waves  from  directions  between  south 
and  west. 

d.  Plan  3,  a  modification  to  Plan  2,  provides  adequate  protection  from 
the  incident  wave  climate  since  die  additional  east  breakwater  exten¬ 
sion  ovo'liqis  the  existing  south  breakwater  and  permits  very  low 
energy  inside  the  harbor.  This  plan  is  also  a  recommended  alterna¬ 
tive. 

e.  Plan  la,  similar  to  Plan  1,  satisfies  the  HQUSACE  design  criteria. 
However,  harbor  p^ormance  rdative  to  the  criteria  may  be  marginal 
in  a  few  bathing  locations,  particularly  along  die  mtisting  wharf. 

Plan  la  is  considered  an  acc^table  alternative.  It  has  an  advantage 
relative  to  Plan  1  in  that  the  south  breakwater  extension  is  farther  west 
and  may  be  less  likely  to  affect  the  Maalaea  Pipeline  surfing  area  East 
of  the  harbor. 

f.  Plan  lb,  with  vertical  sheet-pile  bulkheads  along  the  east  face  of  the 
center  mole,  satisfies  the  critoion  for  entrance  channel  protection  but 
fails  to  provide  adequate  protection  of  some  of  the  more  exposed 
harbor  areas.  Plan  lb  is  not  recommended . 
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Table  1 

Percent  Occurrence  of  Wave  Heights  Versus  Direction 


Wav*  Diraetion.  dag  (from  which  waves  approach) 


Wav*  Haight  ft 

135.0 

157.5 

180.0 

202.5 

225.0 

247.5 

270.0 

Total 

0.00-1.00 

0.00 

0.00 

0.01 

0.02 

0.02 

0.02 

0.00 

0.07 

1.01-2.00 

•  0.03 

0.91 

3.97 

4.32 

1.33 

0.28 

10.84 

2.01-3.00 

•  0.07 

0.86 

5.34 

11.05 

4.16 

2.04 

23.52 

3.01-4.00 

0.00 

0.01 

0.28 

4.26 

11.60 

3.81 

2.24 

22.20 

4.01-5.00 

0.00 

0.00 

0.04 

2.28 

8.63 

2.28 

1.88 

15.11 

5.01-6.00 

0.00 

0.00 

0.02 

1.13 

5.14 

1.52 

0.77 

8.58 

6.01-7.00 

0.00 

0.00 

• 

1.81 

3.38 

0.96 

1.22 

7.37 

7.01-8.00 

0.00 

0.00 

0.09 

0.96 

3.64 

0.83 

0.82 

6.34 

8.01-9.00 

0.00 

0.00 

0.00 

0.19 

3.00 

0.21 

0.66 

4.06 

9.01  -t- 

0.00 

0.00 

0.00 

0.00 

0.00 

1.25 

0.65 

1.90 

TOTAL 

0.01 

0.11 

2.21 

19.96 

50.78 

16.37 

10.56 

100.0 

Table  8 

HARBD  Wave  Amplification  Factors 
Existing  Plan,  Wave  Angle  =  247.5  deg 


av*  Pariod.  aac 


9  Til  M3 


2  10.63 


0746  ro.66 


HARBD  Wave  Amplification  Factors 
Plan  1.  Wave  Angle  =  135.0  deg 


Table  1 1 

HARBD  Wave  Amplification  Factors 
Plan  1,  Wave  Angle  =  157.5  deg 


Table  23 

HARBD  Wave  Amplification  Factors 
Plan  2,  Wave  Angle  =  270.0  deg 


Table  27 

HARBD  Wave  Amplification  Factors 
Plan  3.  Wave  Angie  =  202.5  deg 


Table  33 

HARBO  Wave  Amplification  Factors 
Plan  1a.  Wave  Angle  =  180.0  deg 


Wave  amplitude  is  below  significance  for  tabulation. 


0.01 

0.02 

0.01 

0.02 

HARBD  Wave  Amplification  Factors 
Plan  1a,  Wave  Angle  =  225.0  deg 


Basin  19 


i  fo:94  0.97 


2  0.82  0.81 


0.13  0.18 


0.07  0.11 


0.61  0.61  0.44 


0.46  0.43  0.30 


6 


0.13  10.13  0.10 


0.08  0.08  0.05 


Table  4r 

HARBO  Wave  Amplification  Factors 
Plan  1b.  Wave  Angle  =  180.0  deg 


Wave  Period,  sec 


9 


0.69 


0.42 


0.20 


0.11 


0.08 


IW 


IQ 


i  16 

0.20 

0.26 

ran 

17 

0.12 

0.11 

0.09 

18 

0.05 

0.04 

0.05 

19 

0.03 

0.04 

0.03 

10.07 

0.08 

0.04 

0.07 

0.12 

0.15 

0.05 

0.08 

W 

.31  1 

0.1 1 

Wave  amplitude  is  below  significance  for  tabulation. 


Table  42 

HARBD  Wave  Amplification  Factors 
Plan  1b.  Wave  Angie  =  225.0 


Table  45 

HARBD-SHALWV  Wave  Heights  Exceeding  HQUSACE  Criteria 
(Deepwater  Wave  Conditions) _  _ _ 


Deepwater 

HARBD 

SHALVW 

Period 

Height 

Height 

Amp. 

Height 

Isec) 

(ft) 

(ft» 

factor 

(HI 

Existing  Plan  -  1  ft  Criteria 


Existing  Plan  •  2  ft  Criteria 


15 

2.02 

3.00 

17 

2.07 

3.00 

11 

2.05 

4.00 

17 

2.04 

4.00 

9 

2.01 

4.00 

9 

2.07 

8.00 

*  Daspwatar  wave  heights  hetwsan  1-9  ft  do  not  oxcaod  HQUSACE  criteria  for  this  condi¬ 
tion. 


Table  47 

HARBD'SHALWV  Wave  Heights  Exceeding  HQUSACE  Criteria' 
(Deepwater  Wave  Conditions)  _ 


Oaapwatar 

HARBD 

SHALWV 

Period 

Haight 

Haight 

Amp. 

Haight 

(aae) 

(ftl 

(ftl 

factor 

IftI 

Propoaad  Plan  2  -  1  ft  Critaria 


Proposed  Plan  2  ■  2  ft  Criteria 


1S7.5 

11 

2.09 

4.00 

1.20 

1.74 

1 

180.0 

9 

2.68 

3.00 

1.56 

1.72 

1 

11 

2.35 

3.00 

1.45 

1.62 

1 

13 

2.37 

4.00 

1.16 

2.04 

1 

17 

2.05 

5.00 

0.94 

2.18 

1 

20 

2.04 

6.00 

0.91 

2.24 

1 

202.5 

13 

2.29 

4.00 

1.55 

1.48 

1 

20 

2.07 

5.00 

1.51 

1.37 

1 

225.0 

9 

2.03 

7.00 

0.87 

2.34 

1 

11 

2.04 

7.00 

0.94 

2.17 

1 

Daap water  wave  heights  between  1-9  ft  do  not  exceed  HQUSACE  criteria  for  this  corxfi- 
tion. 


Table  49 

HARBD-SHALWV  Wave  Heights  Exceeding  HQUSACE  Criteria 
(Deepwater  Wave  Conditions)  _ 


Deepwater 

HARBD 

SHALVW 

Basin 

Direction 

Period 

Height 

Height 

Amp. 

Height 

Num. 

Idag) 

leac) 

Ift) 

mi 

factor 

mi 

# 

Proposed  Plan  1  a  -  2  ft  Criteria 


Table  51 

Percent  Occurrence  of  Wave  Heights  Versus  Direction* 
Existing  Plan  -  Wave  Heights  Exceeding  1  ft _ 


Wav*  Diraetion,  d«g  (from  which  wav*s  approach) 


1S7.&  180.0  202.S  22S.0  247.S  270.0 


0.00-1.00 


1.01-2.00 


2.01-3.00 


3.01-4.00 


4.01-5.00 


5.01-6.00 


6.01-7.00 


7.01-8.00 


8.01-9.00 


9.01  + 


TOTAL 


Pflrcent  occurrenca  is  below  sigruficance  for  tabulation. 


1.68 

2.06 

0.34 

0.77 

2.89 

0.36 

0.19 

2.34 

0.13 

1.25  0.65 


5.94  12.00  2.08  0.65 


Table  52 

Percent  Occurrence  of  Wave  Heights  Versus  Direction* 
Existing  Plan  •  Wave  Heights  Exceeding  2  ft 


157.5  180.0  202.5  225.0  247.5  270.0 


Wave  Direction,  dag  (from  which  waves  approach! 

Wav*  Haight  - 

ft  135.0 


0.00-1.00 


1.01-2.00 


2.01-3.00 


3.01-4.00 


4.01-5.00 


5.01-6.00 


6.01-7.00 


7.01-8.00 


8.01-9.00 


9.01 


TOTAL 


Paroant  occurranca  is  balew  aigitificanca  for  tabulation. 


Table  55 

Percent  Occurrence  of  Wave  Heights  Versus  Direction' 
Plan  2  -  Wave  Heights  Exceeding  1  ft 


Wav*  Oiraction,  dag  (from  which  wavas  approach) 


Wave  Haight 
ft 


2.01-3.00 


3.01-4.00 


4.01-5.00 


5.01-6.00 


6.01-7.00 


7.01-8.00 


8.01-9.00 


9.01  + 


TOTAL 


0.32 


0.06  2.09 


0.03  1 .40 


180.0  I  202.5  225.0  247.5  270.0  Total 


0.32 


0.91 

2.79 

1.18 

1.91 

0.70 

2.11 

0.19 

1.97 

0.50  6.47 


2.16 


1.43 


3.70 


3.09 


2.90 


2.16 


1.25  0.65  1.90 


1.25  0.65  17.66 


Parcant  occurrence  is  below  significance  for  tabulation. 


Table  56 

Percent  Occurrence  of  Wave  Heights  Versus  Direction* 
Plan  2  •  Wave  Heights  Exceeding  2  ft 


Wave  Oiraction.  dag  (from  which  wavas  approach) 

Wave  Haight  '  I  "  I - - -  ■  - 

ft  135.0  157.5  180.0  202.5  225.0  247.5  270.0  Total 


2.01-3.00 


3.01-4.00 


4.01-5.00 


6.01-6.00 


6.01-7.00 


7.01-8.00 


8.01-9.00 


9.01  + 


TOTAL 


Parcant  occurranea  is  balow  significance  (or  tabulation. 


0.32 


2.16 


1.29 


0.66 


1.54 


1.69 


1.71 


1.25  0  65  1.90 


1.25  0.65  11.27 


Table  59 

Percent  Occurrence  of  Wave  Heights  Versus  Direction 

Plan  1  a  -  Wave  Heights  Exceeding  1  ft  _ _ 


Wave  Direction,  deg  (from  which  wavea  approach) 

Wava  Haight  |  I  I 

ft  135.0  157.5  180.0  202.5  225.0  247.5  270.0  Total 


3.01 

-4.00 

4.01 

-5.00 

5.01 

-6.00 

6.01 

-7.00 

7.01 

-8.00 

8.01 

-9.00 

9.01  + 


TOTAL 


1.25  0.65  1.90 


0.03  0.97  6.97  1.38  0.65  10.00 


Percent  occurrence  is  below  significance  for  tabulation. 


Table  60 

Percent  Occurrence  of  Wave  Heights  Versus  Direction* 
Plan  la  -  Wave  Heights  Exceeding  2  ft 


Wave  Direction,  deg  (from  which  wave*  approach) 

Wave  Height  i  ,  .  I 

ft  135.0  157.5  180.0  202.5  225.0  247.5  270.0  Total 


2.01-3.00 


3.01-4.00 


4.01-5.00 


5.01-6.00 


6.01-7.00 


7.01-8.00 


8.01-9.00 


9.01  + 


TOTAL 


0.00 


0.00 


0.00 


0.01 


0.35 


1.18 


1.52 


1.25  0.65  1.90 


1.25  0.65  4.96 


Percent  occurrence  is  below  significance  for  tabulation. 


Table  61 

Percent  Occurrence  of  Wave  Heights  Versus  Direction* 
Plan  1  b  -  Wave  Heights  Exceeding  1  ft 


Wav*  Diraction,  dag  (from  which  wavaa  approachl 
IVava  Haight  |  |  III 

ft  13S.0  1S7.S  180.0  202.5  225.0  247.5  270  0  Total 


1.25  0.65  1.90 


0.05  3.96  12.76  1.47  0.65  18.89 


Percent  occurrence  is  below  significance  for  tabulation. 


Table  62 

Percent  Occurrence  of  Wave  Heights  Versus  Direction* 
Plan  1b  •  Wave  Heights  Exceeding  2  ft 


Wave  Oiraetion,  deg  (from  which  waves  approachl 
Wave  Height  — ~~*~- 
ft  135.0 


157.5  180.0  202.5  225.0  247.5  270.0  Total 


2.01-3.00 


3.01-4.00 


4.01-5.00 


5.01-6.00 


6.01-7.00 


0.00  0.16  2.81 


Percent  occurrence  is  below  significance  for  tabulation. 


1.25  0.65  1.90 


1.25  0.65  4.87 


Table  63 

Summary  of  Percent  Occurrence  of  Wave  Heights 


Location 

Percent  of  Time  Criterion  is  Exceeded 

HQUSACE 

Criterion 

Existing 

Plan  1 

Plan  2 

Plan  3 

Plan  la 

Plan  1b 

Berthing  areas 
(1  ft  criterion) 

<  10.0 

21.4 

6.1 

17.7 

B 

10.0 

18.9 

Entrance  Channel 
(2  ft  criterion) 

<  10.0 

9.6 

11.3 

B 

B 
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Table  A-1 

Percent  Occurrence  of  Height  and  Period  by  Direction 
Wave  Direction  =  135.0  deg  (from  which  waves  approach) 


Wave  Height 
ft 

Peak  Period  (seel  |j 

9 

11 

13 

15 

17 

20 

Total 

0.00-1.00 

0.00 

1.01-2.00 

a 

• 

• 

2.01-3.00 

a 

a 

3.01-4.00 

0.00 

4.01-5.00 

0.00 

5.01-6.00 

0.00 

6.01-7.00 

0.00 

7.01-8.00 

0.00 

8.01-9.00 

0.00 

9.01  -h 

0.00 

Total 

0.00 

0.00 

0.00 

0.00 

1 

• 

0.01 

Percent  occurrence  is  below  table  resolution. 


Table  A-2 

Percent  Occurrence  of  Height  and  Period  by  Direction 
Wave  Direction  =  1 57.5  deg  (from  which  waves  approach) 


Wave  Height 
ft 

P«ak  Pariod  laac)  || 

9 

11 

13 

15 

17 

20 

Total 

0.00-1.00 

0.00 

1.01-2.00 

• 

0.01 

0.02 

0.03 

2.01-3.00 

0.03 

0.01 

0.03 

• 

0.07 

3.01-4.00 

0.01 

0.01 

4.01-5.00 

0.00 

5.01-6.00 

0.00 

6.01-7.00 

0.00 

7.01-6.00 

0.00 

8.01-9.00 

0.00 

9.01 

0.00 

Total 

0.00 

0.00 

0.04 

0.01 

0.04 

0.02 

0.11 

*  Percent  occurrence  m  below  table  resolution. 
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Table  A-3 

Percent  Occurrence  of  Height  and  Period  by  Direction 
Wave  Direction  =  180.0  deg  (from  which  waves  approach) 


Wava  Haight 
ft 


Paak  Pariod  laac) 


9 

11 

13 

15 

17 

20 

• 

1 

0.01 

0.01 

0.02 

0.03 

0.14 

0.52 

0.19 

0.09 

0.16 

0.07 

0.06 

0.39 

0.09 

0.03 

0.02 

, 

• 

0.18 

0.04 

0.03 

« 

• 

0.02 

• 

0.09 

g 

0.22 

0.20 

0.13 

0.21 

1.12 

0.33 

0.00-1.00 _ ; _ ; _ ■ 

1.01-2.00 _ 0.01  0.02  0.03 

2.01-3.00 _ 0.09  0.16  0.07 

3.01-4.00 _ 0.03  0.02 

4.01-5.00 _ ; _ ; _ 0.03 

5.01-6.00 _ ; _ ; _ . 

6.01-7.00 _ ^ _ ^ _ . 

7.01-8.00 _ 0.09 _ ■ 

8.01-9.00 _ _ ^ _ . 

9.01  + _ ; _ ; _ . 

Total _ 0.22  0.20  0.13 

Percent  occurrence  is  below  table  resolution 


Table  A-4 

Percent  Occurrence  of  Height  and  Period  by  Direction 
Wave  Direction  =  202.5  deg  (from  which  waves  approach) 


Wave  Height 
ft 


Peak  Pariod  (sec) 


9 

n 

13 

15 

17 

20 

1.18 

0.66 

0.51 

1.70 

0.90 

0.39 

0.47 

0.07 

0.37 

2.09 

1.12 

0.14 

0.15 

0.10 

1.25 

0.61 

0.17 

0.07 

0.18 

0.66 

0.21 

• 

0.08 

1.04 

0.63 

0.06 

0.09 

0.50 

0.26 

0.11 

0.19 

2.32 

0.60 

1.53 

9.25 

4.76 

1.30 

2.01  -3.00  1.18  0.66  0.51 

3.01-4.00  0.47  0.07  0.37 

4.01-5.00 _ 0.15 _ 0.10 

5.01-6.00 _ 0.07 _ 0.18 

6.01-7.00 _ 0.08 _ ; _ 

7.01-8.00 _ 0.09 _ ■ 

8.01-9.00 _ ^ _ ^ _ ■ 

9.01 -f _ ; _ ; _ ■ 

Total _ 2.32  0.60  1.53 

Percent  occurrence  is  below  table  resolution. 


1.13 

1.81 

0.96 

0.19 

0.00 
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Table  A-5 

Percent  Occurrence  of  Height  and  Period  by  Direction 
Wave  Direction  =  225.0  deg  (from  which  waves  approach) 


Peak  Period  (sec) 


W«v»  might 
ft 

9 

11 

13 

15 

17 

20 

Total 

2.01-3.00 

1.46 

2.92 

4.31 

1.46 

0.84 

0.06 

1 1.05 

3.01-4.00 

1.02 

2.37 

5.44 

2.05 

0.68 

0.04 

11.60 

4.01-5.00 

0.90 

1.63 

2.82 

1.98 

1.24 

0.06 

8.63 

5.01-6.00 

0.49 

1.39 

2.30 

0.81 

0.15 

• 

5.14 

6.01-7.00 

0.34 

0.16 

1.41 

1.15 

0.30 

0.02 

3.38 

7.01-8.00 

0.41 

0.68 

1.02 

0.73 

0.73 

0.07 

3.64 

8.01-9.00 

0.88 

0.64 

0.45 

0.54 

0.48 

0.01 

3.00 

9.01  -I- 

0.00 

Total 

6.32 

10.90 

19.38 

9.32 

4.56 

0.30 

50.78 

Percent  occurrence  is  below  table  resolution. 


Table  A-6 

Percent  Occurrence  of  Height  and  Period  by  Direction 
Wave  Direction  =  247.5  deg  {from  which  waves  approach) 

Peak  Period  (cec) 


Wave  Heisht 
ft 

9 

11 

13 

15 

17 

20 

Total 

2.01-3.00 

0.73 

2.22 

0.78 

0.31 

0.08 

0.04 

4.16 

3.01-4.00 

0.64 

1.87 

0.69 

0.51 

0.08 

0.02 

3.81 

4.01-5.00 

0.35 

1.22 

0.38 

0.31 

0.01 

0.01 

2.26 

5.01-6.00 

0.55 

0.68 

0.06 

0.17 

0.05 

0.01 

1.52 

6.01-7.00 

0.34 

0.41 

0.18 

0.03 

• 

0.96 

7.01-8.00 

0.36 

0.28 

0.17 

« 

0.01 

0.83 

8.01-9.00 

0.13 

0.07 

0.01 

0.21 

9.01  + 

1.12 

0.06 

0.07 

• 

1.25 

Tctal 

4.40 

7.52 

2.43 

1.63 

0.29 

0.10 

16.37 

Percent  occurrence  ie  below  table  resolution. 
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Table  A-7 

Percent  Occurrence  of  Height  and  Period  by  Direction 

Wave  Direction  =  270.0  deg  (from  which  waves  approach) 

Wav«  Haight 
ft 

Pttak  Pariod  <«acl  | 

9 

11 

13 

IS 

17 

20 

Total 

2.01-3.00 

0.39 

0.86 

0.60 

0.16 

0.02 

■ 

2.04 

3.01-4.00 

0.47 

0.88 

0.80 

0.04 

0.03 

0.02 

2.24 

4.01-5.00 

0.11 

1.03 

0.61 

0.12 

• 

« 

1.88 

5.01-6.00 

0.31 

0.32 

0.11 

0.03 

■ 

0.77 

6.01-7.00 

0.07 

0.53 

0.51 

0.10 

a 

• 

1.22 

7.01-8.00 

0.17 

0.19 

0.44 

0.02 

0.82 

8.01-9.00 

0.12 

0.54 

0.66 

9.01  + 

0.22 

0.32 

.09 

.02 

0.65 

Total 

1.66 

4.72 

3.48 

0.57 

0.09 

0.04 

10.56 

Percent  occurrence  is  below  table  resolution. 
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Appendix  B 
Notation 


a  Wave  amplitude  function 
ao  Incident  wave  amplitude 
c  Wave  celerity 

Cg  Group  celerity 

g  Gravitational  acceleration 
H  Wave  height 

h  Water  depth 

i  imaginary  unit  =  (-1)*^ 

K,  Reflection  coefficient 

K  Wave  number  =  2ir/L 

L  Wavelength 

n  Independent  variable  in  the  direction  of  the  unit  vector 
s  Wave  phase  function 
T  Wave  period 

X  Horizontal  coordinate 

y  Horizontal  coordinate 

a  Reflective  component  of  absorbing  boundary 

Dimensionless  bottom  friction  coefficient 
7  Phase  difference  between  bottom  friction  and  flow  velocity 
6  Wave  approach  angle 
X  Complex  bottom  friction  factor 
IT  3.14159 . 

2)  Radian  wave  frequency,  intrinsic  wave  frequency 
<t>  Velocity  potential 
d  Partial  differentiation 

V  Gradient  operator  in  two  dimensions  =  (d/3x  +  didy) 
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